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Abstract. We provide two certificates of convexity for arbitrary basic closed 
semi-algebraic sets of R". The first one is based on a necessary and suffi- 
cient condition whereas the second one is based on a sufficient (but simpler) 
condition only. Both certificates are obtained from any feasible solution of a 
related semidefinite program and so in principle, can be obtained numerically 
(however, up to machine precision). 



1. Introduction 

With R[x] being the ring of real polynomials in the variables x\, . . . , x n , consider 
the basic closed semi-algebraic set K C W 1 defined by: 

(1.1) K:= {x&W 1 : 9j (x) > 0, j = l,...,m} 

for some given polynomials gj el[i],j = l,...,m. 
By definition, K C R™ is convex if and only if 

(1.2) i,!/eK => Xx + (l-X)yeK VAe[0,l]. 

The above geometric condition does not depend on the representation of K but 
requires uncountably many tests and so cannot be checked in general. 

Of course concavity of gj for every j = 1, . . . , m, provides a certificate of convex- 
ity for K but not every convex set K in (jl.ll) is defined by concave polynomials. 
Hence an important issue is to analyze whether there exists a necessary and suf- 
ficient condition of convexity in terms of the representation (jl.ip of K because 
after all, very often (jl.ip is the only information available about K. Moreover, a 
highly desirable feature would be that such a condition can be checked, at least 
numerically. 

In a recent work [3] , the author has provided an algorithm to obtain a numerical 
certificate of convexity for K in (|1.1[) by using the condition: 

(1.3) (V gj {y),x- y) > 0, Vx,y G K with 9j (y) = 0, 

which is equivalent to (|1.2[) provided that SlateiQ condition holds and the nondegen- 
eracy condition Vgj(y) ^ holds whenever y £ K and gj(y) = 0. This certificate 
consists of an integer pj and two polynomials 6{, 9% G M[x, y] for each j = 1, . . . , m, 
and their characterization obviously implies that (jl.3l) holds true and so K is convex 
(whence the name certificate); see Lasserre [H Corollary 4.4]. More precisely, for 
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every j = 1, . . . , m, define the 2m + 1 polynomials hi G R[x, y] by ft^(a;, y) = gi{x), 
h m+ e(x,y) = ge{y) for every I = l,...,m, and ft 2 m+i(a;, y) = -gj{y). The pre- 
ordering Pj C M.[a;,y] generated by the polynomials (/i^) C ?/] is defined by: 

(j j G E[x,y] 

I JC{i,...,2m+i} \ieJ / 

where y] C M[x, y] is the set of polynomials that are sums of squares (in short 
s.o.s.), and where by convention, Yiiej he(x, y) = 1 when J = 0. Then by a direct 
application of Stcngle's Positivstellensatz [1, Theor. 4.4.2, p. 92] (more precisely, 
a Nichtnegativstellensatz version) (|1.3[) holds if and only if 

(1.5) e{(x,y)(W 9j (y),x-y) = ({V 9j (y),x - y)f^ + 4{x,y), 

for some integer pj and some polynomials O{,02 G Pj- In addition, bounds (p, d) 
are available for the integer pj and the degrees of the s.o.s. polynomials <jj ap- 
pearing in the definition (|1.4p of polynomials 0{ , 9 3 2 G Pj , respectively. Observe 
that in (|1.5[) one may replace j?^ with the fixed bound p (multiply each side with 
((V gj(y),x — y)) 2 ( p ~ Pi ') and take d := d + p. Next, recall that s.o.s. polynomials 
of bounded degree can be obtained from feasible solutions of an appropriate semi- 
definite progran^ (see e.g. [3]). Hence, in principle, checking whether (|1.5|) has a 
feasible solution (#^,#2) reduces to checking whether a single semidefmite program 
has a feasible solution. 

And so, when both Slater and the nondegeneracy condition hold, checking whether 
K is convex reduces to checking if each of the semidefmite programs associated with 
(jl.5j) . j = 1, . . . , m, has a feasible solution. When K is convex, the 2m polynomials 
91,02 £ ^[^Z/L J — 1, • • - ,"2, provide the desired certificate of convexity through 
(ll.3|) ; see [U Corollary 4.4]. However, it is only a numerical certificate because it 
comes from the output of a numerical algorithm, and so subject to anavoidable nu- 
merical inaccuracies. Moreover, the size of each semidefmite program equivalent to 
(|1.5j) is out of reach for practical computation, and in practice, one will solve a semi- 
definite program associated with ()1.5[) but for reasonable bounds (p' 7 d') (p,d), 
hoping to obtain a solution when K is convex. An alternative and more tractable 
certificate of convexity using quadratic modules rather than preorderings is also pro- 
vided in [4j Assumption 4.6], but it only provides a sufficient condition of convexity 
(almost necessary when K is compact and satisfies some technical condition). 

The present contribution is to provide a certificate of convexity for arbitrary 
basic closed semi-algebraic sets i.e., with no assumption on K. This time, 

by certificate we mean an obvious guarantee that the geometric condition (|1.2[) 
holds true (instead of (|1.3p in 4 ). To the best of our knowledge, and despite 
the result is almost straightforward, it is the first of this type for arbitrary basic 
closed semi-algebraic sets. As in [4] our certificate also consists of two polynomials 
of M.[x, y) and is also based on the powerful Stengle's Positivstellensatz in real 
algebraic geometry. In addition, a numerical certificate can also be obtained as the 
output of a semidefinite program (hence valid only up to machine precision). We 



1 A semidefinite program is a convex optimization problem with the nice property that it 
can be solved efficiently. More precisely, up to arbitrary fixed precision, it can be solved in time 
polynomial in its input size. For more details on semidefinite programming and its applications, 
the interested reader is referred to e.g. [H]. 
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also provide another certificate based on a simpler characterization which now uses 
only a sufficient condition for a polynomial to be nonnegative on K; so in this case, 
even if K is convex, there is no guarantee to obtain the required certificate. Finally, 
we also provide a sufficient condition that permits to obtain a numerical certificate 
of non convexity of K in the form of points x, y G K which violate (|f ,2[) . 

2. Main result 
Observe that in fact, (|f ,2p is equivalent to the simpler condition 

(2.1) x,yeK => (x + y)/2eK. 

Indeed if K is convex then of course (|2.1[) holds. Conversely, if K is not convex 
then there exists x, y G K and < A < 1 such that z := x + X(y — x) is not in K. As 
K is closed, moving on the line segment [x, y] from x to y, there necessarily exist 
x € K (the first exit point of K) and y G K (the first re-entry point in K), with 
x 7^ y. Thus, as x and y are the only points of [x, y] contained in K, the mid-point 
z := (x + y)/2 is not contained in K. 

Given the basic closed semi- algebraic set K defined in (jl.ip . let K :— K x K C 
R™ x E™ be the associated basic closed semi-algebraic set defined by: 

(2.2) K := {(x,y) : y) > 0, j = 1, . . . , 2m }, 
where: 

(2.3) (x,y) i-> gj(x,y) := g-j(x), j = l,...,m 

(2.4) (x,y) t-^ gj(x,y) := g m -j(y), 3 = m + 1, . . . , 2m, 

and let P(g) C R[x, y] be the preordering associated with the polynomials (gj) that 
define K in lpT2|) . i.e., 

(2-5) P(s) := | ^ ^(n^) : 0jeE[aj,y] I, 

[ JC{L...,2m} VfeSJ / J 

where ~S[x, y] C R[x, y] is the set of s.o.s. polynomials. Our necessary and sufficient 
condition of convexity is a follows. 

Theorem 2.1. Let K C K™ fee £/ie 6asic closed semi- algebraic set defined in 
Then K is convex if and only if for every j — l,...,m, there exist polynomials 
o~j, hj G P(g) and an integer pj G N such that: 

(2.6) a j (x,y)g j ((x + y)/2) = 9j ((x + y)/2) 2 ^ + hj(x, y), Vx,yGM". 

Proof. The set K is convex if and only if (|2.1|) holds, that is, if and only if for every 
j = 1, ...,m, 

(2.7) gj ({x + y)/2) >0 V(x,y)GK. 

But then (|2.6[) follows from a direct application of Stengle's Positivstellensatz [U 
Theor. 4.4.2, p. 92] to (|2.7j) (in fact, a Nichtnegativstellensatz version). □ 

The polynomials o~j,hj G -P(s), j = 1, • ■ ■ ,m, obtained in (|2.6I) indeed provide 
an obvious certificate of convexity for K. This is because if (|2.6[) holds then for 
every x, y G K one has o~j{x, y) > and hj(x, y) > because <Xj, /ij G P{g)] and so 
(jj (x, y)gj((x + y) /2) > 0. Therefore if Oj (x, y) > then gj ((x + y)/2) > whereas 
if o-j(x, y) — then ^ ((x + y)/2) 2pj = which in turn implies gj((x + y)/2) = 0. 
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Hence for every j = 1, . . . , m, gj((x + y)/2) > for every x, y G K, that is, (|2.ip 
holds and so K is convex. 

A numerical certificate of convexity. Again, as (|2.6p is coming from Stengle's 
Positivstellensatz, bounds (p, d) are available for the integer pj and the degrees of 
the s.o.s. polynomials <pj appearing in the definition (|2.5p of polynomials <jj, hj G 
P(g), respectively. Hence, with same arguments as in the discussion just after (j 1 . 5[) . 
checking whether (|2.6|) holds reduces to check whether some (single) approriately 
defined semidefinite program has a feasible solution. 

Hence checking convexity of the basic closed semi-algebraic set K reduces to 
checking whether each semidefinite program associated with (|2.6I) . j = 1, . . . ,m, 
has a feasible solution, and any feasible solution (Xj, hj e P(g) of (|2.6| . j = 1, . . . , m, 
provides a certificate of convexity for K. However the certificate is only "numerical" 
as the coefficients of the polynomials o~j, hj are obtained numerically and are subject 
to anavoidable numerical inaccuracies. Moreover, the bounds (p, d) being out of 
reach, in practice one will solve a semidefinite program associated with (|2.6[) but 
for reasonable bounds (p' ,d') <C (p,d), hoping to obtain a solution when K is 
convex. 

2.1. An easier sufficient condition for convexity. While Theorem l2.1l provides 

a necessary and sufficient condition for convexity, it is very expensive to check 
because for each j = 1, . . . , m, the certificate of convexity Uj, hj E P(g) in (|2.6p 
involves computing 2 x 2 2m = 2 2m+1 s.o.s. polynomials <f>j in the definition (|2.5p 
of <ij and hj . However, one also has the following sufficient condition: 

Theorem 2.2. Let K C K™ he the basic semi- algebraic set defined in Then 
K is convex if for every j = 1, . . . , m: 

m 

(2.8) gj ((x + y)/2) = o- (x, y) + J2 a k( x ' v)9k(x) + ^{{x, y) g k (y), Va, y G R n , 

fc=i 

for some s.o.s. polynomials G E[x, y]. 

Proof. Observe that if (|2.8p holds then gj((x + y)/2) > for every j — 1, . . . , m 
and all (x,y) G K; and so K is convex because (|2.ip holds. □ 

Again, checking whether (|2.8I) holds with an apriori bound 2d on the degrees of 
the s.o.s. polynomials o~ 3 k ,il) 3 k , reduces to solving a semidefinite program. But it 
now only involves 2m + 1 unknown s.o.s. polynomials (to be compared with 2 2m+1 
previously). On the other hand, Theorem 12.21 only provides a sufficient condition, 
that is, even if K is convex it may happen that (|2.8p does not hold. 

However, when K is compact, convex, and if for some M > the quadratic 
polynomial x h-> M — \\x\\ 2 can be written 

m 

(2.9) M-\\x\\ 2 = <j (x) + J2°k(x)gj(x), 

fe=i 

for some s.o.s. polynomials (tTfc) C S[x], then (|2.8|1 is almost necessary because for 
every e > 0: 

m 

(2.10) 9j {(x + y)/2) + e = 4 e {x, y) + £ o^Jx, y)g k (x) + ^ ke (x, y)g k (y), 

fe=i 
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for some s.o.s. polynomials cr^,"^ 6 S[x, j/]. Indeed, consider the quadratic 
polynomial 

(x,y)^A(x,y) := 2M - \\x\\ 2 - \\y\\ 2 . 
From p.9p . A belongs to the quadratic module Q(g) C M.[x, y] generated by the 
polynomials g k that define K, that is, the set 

Q{g) ■= |CTo(a;, y) + Y^a k (x 1 y)g k (x) + ip k (x,y) g k (y) : a k , ip k e E[z, y] j . 

In addition, its level set {(x, y) : A(x, y) > 0} is compact, which implies that Q{g) 
is Archimedean (see e.g. [7]). Therefore, as gj((x + y)/2) + e > on K, (|2.10|) 
follows from Putinar's Positivstellensatz [6]. 

2.2. A certificate on non-convexity. In this final section we provide a numerical 
certificate of non convexity of K when the optimal value of a certain scmidcfinite 
program is strictly negative and some moment matrix associated with an optimal 
solution satisfies a certain rank condition. 

Given a sequence z = {z a p) indexed in the canonical basis (x a y@) of M.[x, y], let 
L z : M.[x, y] — > R be the linear functional: 

a, p a,P 

and as in [3], the moment matrix M s (z) associated with z is the real symmetric 
matrix with rows and columns indexed in the the canonical basis (x a y^) and with 
entries 

M s (z)((a,P),(a',0')) = z ( a + a >)(/3+i3>) 

for every (a, /3), (a',/3') 6 Nf n , where N™ := {a 6 N n : "i < »}< 

Similarly, with a polynomial (x,y) i-* 0(x,y) = J2 a p x a y^ , the localizing 
matrix M s (8z,) associated with and z, is the real symmetric matrix with rows 
and columns indexed in the the canonical basis (x a yP) and with entries 

M s (6z)((a,P),(a',P')) = 9 a »p» z^ a+ct , +a n^ 0+0 , + ^, 

a",0" 

for every (a, /3), (a',0') 6 N^ n . 

Let «fe := [(deg <?fc)/2] , k = 1, . . . , 2m, and for every j = 1, . . . , m, and s > v := 
maxfc v k , consider the scmidcfinite program: 

p js = mm J L z ( i g J ((x + y)/2)) 

s.t. M s (z) b 

M s ^ Vk (g k z) b0, fc = l,...,2m 
z = 1, 



(2.11) 



where for a real symmetric matrix A, the notation A y stands for A is positive 
scmidcfinite. The semidefinite program (|2.11[) is a convex relaxation of the global 
optimization problem 

g* := min { 9j ((x + y)/2) : (x,y) € K } 

and so p J5 < .g* for every 5 > v. Moreover, pj S t 5* as s -— > 00; for more details see 
e.g. 0. 
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Theorem 2.3. Let K C K n be as in (1. 1\) and let z be an optimal solution of the 
semidefinite program h2.11\) with optimal value pj s . If pj s < and 

(2.12) rank M s (z) = rank M a _„(z) (=: t) 

then the set K is not convex and one may extract t points (x(i),y(i)) £ K, i = 
1, . . . , t, such that 

9j {{x{i) + y{i))/2) < 0, Vi = l,...,t. 

Hence each mid-point (x(i) + y(i))/2 ^ K is a certificate that K is not convex. 

Proof. By the flat extension theorem of Curto and Fialkow [2] (see also Laurent 
[5]), the rank condition (|2.12[) ensures that z is the moment sequence of a t-atomic 
probability measure fj, supported on K. That is: 

z af > = I x a y$d^ V(a,j8) eN£. 
Jk 

Let (x{i),y{i))\ =1 C K be the support of /i which is a positive linear combination 
of Dirac measures &(x{i),y{i)) with positive weights (7^) such that X^7« = 1- Then 

.g* > P JS = £-(ffj((a: + y)/2)) = / ff j((a; + y)/2) d M 

t 

= Z)7i5j((«W+l/W)/2) 
»=i 

i=l 

which shows that p J;s = g* and so, gj((x(i) + y(i))/2) = g* for every i — l,...,t. 
But then the result follows from pj s < 0. □ 
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